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Abstract 


We analyze hyperon and nucleon mass shifts in nuclear matter using chi¬ 
ral perturbation theory. Expressions for the mass shifts that include strong 
interaction effects at leading order in the density are derived. Corrections 
to our results are suppressed by powers of the Fermi momentum divided by 
either the chiral symmetry breaking scale or the nucleon mass. Our work is 
relevant for neutron stars and for large hypernuclei. 
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In nuclear matter the masses of hyperons and nucleons are shifted from their vacuum 
values by strong interaction effects. Understanding these mass shifts is important for neu¬ 
tron stars ra and for large hypernuclei (see for example ||||). Most approaches to this 
problem rely on the use of phenomenological models 1[I|.|2|.|4| pl.|ti| [Hi].? Jl5|| since the relevant 
strong interaction effects involve QCD in the nonperturbative regime. In this letter we will 
consider mass shifts for the baryon octet using a model independent approach to baryon- 
baryon interactions based on chiral perturbation theory. The price to be paid for model 
independence is that an expansion in momentum is necessary and that our expressions for 
some of the mass shifts involve parameters in the heavy baryon chiral Lagrangian that have 
not yet been determined from experiment. However, there are mass relations that hold 
independent of these parameters in analogy with the Gell-Mann-Okubo mass formula for 
SU(3) breaking. 

We study changes in the masses of the baryon octet (n,p, A, E, S) that arise from their 
strong interactions with a degenerate Fermi gas of neutrons at density S 71 ' 1 = (p^) 3 / 37T 2 and 
protons at density S p ' ) = (p b ' > ) 3 /3tt 2 (i.e., nuclear matter). These interactions are determined 
by a heavy-baryon chiral Lagrangian that is consistent with the SU{3) B x SU(3) b chiral 
symmetry of QCD. An expansion in derivatives is appropriate for Fermi momenta pp and 
pp that are small compared with the chiral symmetry breaking scale and the baryon masses. 

In the rest frame of background nuclear matter the energy of a baryon and its three- 
momentum are related (for \p\ « M B ) by E NM = M B + A M B + \p\ 2 /2M B + p B \p\ 2 + ... 
which is to be compared to the free space relation E vac = M B + \p\ 2 /2M B +... where the dots 
denote terms higher order in the momentum expansion. We can combine the effects of the 
medium with the free space mass, M B , by defining M B and M B , with M B = M B + AM B and 
M b = M b ( 1 + 2r) B M B )~ 1 = M b + A Mg allowing us to write E NM = M B + \p\ 2 /2M B + .... 
The quantities AM B and r/ B arise from interactions with the background medium and in 
general 7 ^ M B . We compute the lowest order contribution to the momentum independent 
mass shift A M B and the coefficient of the |p | 2 term r] B for the octet baryons using chiral 
perturbation theory. 

The pseudo-Goldstone boson fields ( 7 r, K, rj) can be written as a 3 x 3 special unitary 
matrix 


E 


2iU 



( 1 ) 


where 


n 


‘7T°/\/2 + p/VQ 7T + 

7 — 7 T 0 /\/2 + r)/\/E 
K- K° 



( 2 ) 


Under chiral SU(3)l x SU(3) b symmetry, E —>■ LT,W where LeSU(3)L and ReSU(3) B . At 
leading order in chiral perturbation theory / can be identified with the pion decay constant 
(/tt ~ 132MeU). When describing the interactions of the pseudo-Goldstone bosons with 
other fields it is convenient to introduce 


£ = expy = a/E 


(3) 
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Under chiral SU(3)l x SU(3)r symmetry 

£ - L= U£1$ 


( 4 ) 


where in general U is a complicated function of L, R and the meson fields II. For transfor¬ 
mations V — L — R in the unbroken SU(3)v subgroup U = V. 

The baryon helds are introduced as a 3 x 3 octet matrix 


B 


■E%/2 + A/V6 £ + p - 

E- -E°/a/2 +A/a /6 n 

E~ S° - a /2a 

L V 


(5) 


that transforms under chiral 577(3 )l x SU(3)r as 


5 -»• UBU^ 


( 6 ) 


We construct a chiral Lagrangian by treating the baryons as heavy static 2-component 
fields. The chiral Lagrangian is written as 

£ = £ (1) + £ (2) + ... , (7) 

where £ (j ) contains 2 j baryon fields. 

£7) is the familiar heavy baryon chiral Lagrangian that gives the interactions of the 
baryon octet with the pseudo-Goldstone bosons. At leading order in chiral perturbation 
theory 

£ 7 ) = TrBjidoBj + iTrB}[V 0l Bj\ 

-DTrB]a jk {A,B k } - FTrB]a jk [A,B k \ , (8) 

with the repeated spin indices j and k summed over 1,2 and the vector and axial-vector 
chiral fields are 

u = +sy?') . (9) 

A. = - « 9 X f ) • (io) 

Nuclear (3 decay and semileptonic hyperon decay give F ~ 0.44 and D ~ 0.81 at tree-level. 
At higher order in chiral perturbation theory there are corrections to Eq. (j||) coming from 
terms with more derivatives and terms with insertions of the light quark mass matrix. 

Interactions between baryons mediated by pseudo-Goldstone boson exchange give rise to 
hyperon mass shifts. At the same order in chiral perturbation theory (i.e., order pp) terms 
in £ with four baryon fields (and no derivatives) also play a role. They are given by 

£ |2) = - - ^IXfljBjBtBj) 

- ^|rr(B»B,)Tr(BjB,) - ^Tr(BtB 3 -)Tr(Bt Bj ) . (11) 
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Factors of 1 /f 2 appear in Eq. (|lTD so that the coefficients Cj are dimensionless. At higher 
order in chiral perturbation theory there are corrections to Eq. (|TT|) coming from terms with 
derivatives and terms with insertions of the light quark mass matrix. There are also contribu¬ 
tions to the baryon mass shifts from operators in C (j> , :') > 2 , involving more baryon fields, 
e.g. Tr(B}BiBjBjB' k Bk), but these are suppressed by additional powers of S n,p ^/(4:irf 3 ) 
over the contribution from operators appearing in Eq. (|ll|) . Note there is no term involving 
Tr(BjBj) Tr(BiBj) in Eq. ([TT|) as it can be eliminated using a Cayley-Hamilton operator 
identity: 

-Tr(BjBjBiBj) + Tr(B|BjB,B i ) - irr(BjBjBjBi) + ^Tr(BtB,B}B s ) 

= ^Tr(B}Bi)Tr(BjBj) - irr(BjBj)Tr(BjBO - irrfBjBjiTrfB.B,) , (12) 


It is easy to understand from spin-flavour considerations why there are only six independent 
four-baryon operators. The SU(3)y decompostion of the product of two baryon octets is 
8 ® 8 = 27 © 10 © 10 © 8 a © 85 © 1 of which the 27, 85 ,1 are symmetric and the 10,10, 8 a 
are antisymmetric under baryon interchange. Fermi statistics require that the baryons in 
the 27, 85 ,1 have spin 5 = 0 and those in the 10,10, 8 a have S = 1 as the operators under 
consideration give rise to S-wave interactions only. Hence, there are only 6 independent 
four-baryon operators that exist at this order in chiral perturbation theory. It is these 
six operators (and pseudo-Goldstone boson exchange) that give the leading contribution to 
S-wave baryon-baryon scattering, such as AN —> AN, in chiral perturbation theory. An 
SU(3)y analysis of hyper on-nucleon scattering can be found in [ pfj . 

The four-baryon operators in Eq. (|TT|) include effects from the repulsive core of the 
baryon-baryon potential. For example, a repulsive spherical well of infinite height and width 
r 0 gives rise to coefficients, Cj, of order AkPtq/Mb- 

Using Eq. (] 8 |) and Eq. (FLj) we End the following one loop expressions for the baryon 
octet mass shifts in nuclear matter: 


A M n 


A M p 


A M-£+ 


AM S - 



d® 



d® 


P 

d ( n > f 

7m 

d^ 

+ JT 

d ( n > ( 

7M 


(F + .D) j (m^/pP) + ^—J^Lx^/pp^) + ci - c 2 + c 5 - c 6 | 

- 2 — /Pf' > ) + 2 Cj + c 2 + 2 C 5 + c 6 | 

-- 2 —/Pj7) + 2 ci + c 2 + 2 c 5 + C 6 j 

| (F + £>) j {rn^/p'f) + ^—P)_x(m v /p { p) + c x - c 2 + c 5 - c 6 | 

— + C 4 + 2 C 5 + Cg j 

j — , ? F - l(m K /p { p) - Cl - 2 c 2 + 2c 5 + c 6 | 

— l(m K /pP) - ci - 2 c 2 + 2c 5 + c 6 | 
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d [p) r c 3 

1 7T + C 4 + 2 C 5 + Cg 


f 2 12 

<^(p) r 1 

Atf=« = { c 3 + 2c 4 + 2c 5 + eg} H —— ^ —c 3 + c 4 + 2 c 5 + eg 

am s - /2 l2 ~o ■ - ■ — ■ -j ' f 
For the A and £° there is a 2 x 2 mass matrix with entries 

d™ { 1 


P 1 ° * ' d p 12 

d H fl „ 1 d (p) r 

c 3 + c 4 + 2c5 + Cg )• H —— |c 3 + 2 c 4 + 2c 5 + Cg} 


( 13 ) 


AM aa = 


{ 1 7 1 13 7 

— (H + 3 F) 2 X{m,K/p^) + 77 c i + x c 2 + T 77 C 3 + -c 4 + 2 c 5 + c 6 
12 o 3 12 o 

{ t 7 r(-D + 3 F) 2 T(m K /pp' 1 ) + -c 4 + -c 2 + —c 3 + -c 4 + 2 c 5 + c 
(12 o 3 12 o 


/\M-£ °s° — 


/ 2 

d® [ 1 

+ T 2 

d( n ) ((D-F ) 2 


/ 2 l 

rf(p) [ (£) - F ) 2 

+ 7 5 l 4 


^ ITT-k/P f'*) — 2 Cl — C 2 + ^ c 3 + 2 c 4 + 2 c 5 + Cg 


(p) \ 1 1 1 

Z{m K /p F ) - -ci - c 2 + — c 3 + — c 4 + 2 c 5 + c 6 


A ~ A d (n) f (D + 3F)(D*F)_ . (n) , ci c 2 5 

A - - EA “ T 5 "! 1X5 * + 2VS + vS " 4V 


4V3 C3 vT 4 , 


d<*> ( (D + 3F)(F - D) , m _ 

/ 2 \ 4V3 ( ,i/Pf > 2i/3 V3 4\/3 A ' V 2 ! 


Cl C 2 5 1 

+ -—7=c 3 H—4=c 4 


In Eq. (|T3| ) and Eq. (|14|) the function Z(x) is defined by 

Z(x) = 1 — 3x 2 + 3x 3 Arcta7i(l/x ) . 


(14) 


(15) 


Naively, corrections to Eq. (|13]) and Eq. (14) are suppressed by powers of the Fermi mo¬ 
mentum divided by the chiral symmetry breaking scale (or the nucleon mass) and by powers 
of the light quark masses divided by the chiral symmetry breaking scale (or the nucleon 
mass). This is essentially a consequence of dimensional analysis as Feynman diagrams with 
more than one loop constructed from the vertices of the Lagrange densities in Eq. (^]) and 
Eq. O contain more powers of 1// which must be compensated by additional factors of p F 
or light quark masses. Corrections to the Lagrange densities in Eq. (^) and Eq. ( |TTD from op¬ 
erators with more derivatives (or insertions of the light quark mass matrix) have coefficients 
with more powers of 1// which again must be compensated by additional factors of p F or 
light quark masses. However, this power counting is not quite correct. Feynman diagrams 
with more baryon propagators than loops have infrared divergent loop integrations over en¬ 
ergy. Furthermore, insertions of the kinetic energy, Cke = TrBjV 2 Bi/2M B , increase the 
number of baryon propagators but not the number of loops (this correction to the Lagrange 
density in Eq. (|8]) cannot be treated as perturbation). The cure for these problems is to sum 
insertions of the kinetic energy changing the baryon propagator to 


- |p|) 


+ 


©GpI -pf) 


p° — p 2 /2M B — ie p°—p 2 /2M B + ie 


(16) 
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Then infrared divergent p° integrations are cut off by the kinetic energy which produces 
factors of M B in the numerator. At /-loops one finds corrections suppressed by a factor of 
order (M B p f /(4vr/ 2 )) £_1 . As a consequence, the leading correction to our results arise from 
two loop diagrams and are suppressed by only rv./ M b Pf /( 47 t f 2 ) (instead of rv./ Pf/ (d^/) 2 )- At 
nuclear density (jpp ~ 350 MeV/c) this factor is not small. Therefore the leading order results 
derived here may receive significant corrections from higher orders in chiral perturbation 
theory. 

We adopt power counting rules where M B is considered to be of the same order as the 
chiral symmetry breaking scale ~ Arcf and two factors of p B are considered to be of the 
same order as a light quark mass (Recall that the squares of the pseudo-Goldstone boson 
masses are proportional to light quark masses). Consequently in Eq. (|13|) and Eq. (|T|) the 
full dependence of the baryon mass shifts on ( Pf/P ) where p = or m v is kept. 

The baryon mass shifts given in Eq. (|13]) and Eq. (|14|) are our principal results. If the 
nuclear matter is not isospin symmetric (i.e., S n> 7 ^ dS p ^) there is A — E° mixing (Eq. ([ID- 
More importantly, the eight octet masses are expressed in terms of six coefficients Cj so there 
are mass relations independent of the coefficients of the four baryon operators. We find that 

2AA/y; 0 S 0 — Ail7y;+ — AAfy;- = 0 , (17) 


which is independent of the composition of the nuclear matter, it is true for an arbitrary 
proton to neutron ratio. A second mass relation is given in terms of the mesonic contribution 

d (n) (A M s+ - A M s _) + c/ (p) (AM h o - AM s _) = ( ° ~P~ (l(m K /pf) - l(m K /pP)) 

(18) 


We see that in the extreme case of neutron matter ( d^ = 0) this relation becomes 
AMs+ = AM=-. There is also a relation involving the E°A mixing term 


f 2 d^ (—3AM aa + 2A M p - 2v / 3AM as - 2AM S - + 3AM s o s o) 
fd (p) (3 AM A a - 2A M n - 2V3 AM A x - 2AM S - + AM s o s o) 

— (D + F) 2 d in)2 + d {p)2 - p n) d ip) (l{m n /pP) - XirnJpP)) 

— (D + F)(D~ 3 F) d (n) d (p) ( l(m K /pP) - l{m K /pf' 

~{D + F) 2 (d^ 2 Z{m K /pP) - d^ 2 X{m K /pf)) 

~\{D- 3 F) 2 d^d^ (l(mJpP) - l(m v /pP)) 


(19) 


When dW = ci- n ) the background nuclear matter is isospin symmetric leading to baryons in 
the same isonmltiplet receiving the same mass shift and E°A mixing vanishing. 

Weinberg pT7| , p~8| has extracted values for c 2 + c 6 = / 2 Ct and 2ci + c 2 + 2c 5 + c 6 = 
f 2 Cs from the measured S-wave nucleon scattering lengths . He found Ct ~ (l/45MeV ) 2 
and C B ~ —(l/ 88 MeV ) 2 which give very large nucleon mass shifts at nuclear density ||19|| . 
If Eq. (p~3[) and Eq. (|R|) are to be relevant at nuclear density the coefficients Cs,t must 
be smaller than these values. It may be possible to extract other combinations of the 
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coefficients c* by measuring the S-wave scattering lengths for hyperon nucleon scattering. 
For hyperon-proton strong interactions the relevant spin singlet as=o and spin triplet 05=1 
elastic scattering lengths are 


|as = o(Ap -*• A p) 
\a s= i(Ap A p) 

|a 5=0 (S + p -► E + p) 
|a s= i(E + p -> E + p) 
\a s =o{^~P S ~p ) 

|a s=1 (E~p -► E _ p) 
|a s=0 (S”p -»• E~p) 
|a s= i(5 _ p -*■ S"p) 

|a s=0 (S°p - E°p) 


P 


0„ . -rO 


P) 


1 

M P M A 

2vr P 

M p + Ma 

1 

m p m a 

2vr/ 2 

Mp + M a 

1 

MpM-% 

2vr/ 2 

Mp + My, 

1 

M p My 

2vr/ 2 

M p + My 

1 

M p My 

27r/ 2 

Mp + My 

1 

M p My 

27r/ 2 

Mp + My 

1 

MpM : s 

27r/ 2 

M p + M H 

1 

MpM= 

27r/ 2 

M p + Ms 

1 

MpMs 

27r/ 2 

Mp + Ms 

1 

MpMs 

27T/ 2 

M p + Me 


,5 5 1 1 

X C 1 — o C 2 — 7C3 + -7C4 + 2C5 — 2 cg 

0000 
3 3 

| Ci + C2 + —C3 + —C4 + 2 c 5 + 2c6 I 

21 Ci — C2 + C5 — Cq I 
21 — Ci — C 2 + C 5 + eg | 

| — C3 + C4 + 2C5 — 2c 6 | 

| C3 + C4 + 2 c 5 + 2 cq \ 

2 1 — C3 + C4 + C5 — C6 1 

2|c 3 + c 4 + c 5 + c 6 | 

| c 3 — C4 — 2C5 + 2 c 6 | 

IC3 + C4 + 2 c 5 + 2 cq \ 


( 20 ) 


Because of the derivative coupling in Eq. (| 8 |) tree-level pseudo-Goldstone boson exchange 
doesn’t contribute to the scattering lengths. 

The baryon mass shifts given in Eq. (|13D and Eq. (0) may be of use for large hypernuclei 
and for neutron stars. Assuming a large nucleus can be modeled by a degenerate Fermi gas 
replacing one of its neutrons with a A results in a hypernucleus with a mass shift that is 
partly determined by AA4 a — AM n . Study of the S-wave hyperon-nucleon scattering lengths 
may eventually lead to a determination of the coefficients Cj which can then be used to 
predict hyperon mass shifts in neutron stars (where S- p > -c d (n - ) ) and in large hypernuclei. 

The mass that appears in the kinetic energy Mg differs from the vacuum mass Mb 
through momentum dependent interactions with the background nuclear matter that give 
a non-zero value for the coefficients ijb■ At one loop the four-baryon operators in Eq. © 
do not contribute to 7]b as they are independent of the baryon momentum. However, the 
meson graphs do give a calculable contribution to t)b and we fold that 


J. n ) /I 1 V 

>f» = ( 5 (B + F) 2 D (+ -(D - 3 FfD(pP/m„)J 

(p) 

+ F) 2 D(j>P/m„) 

Vp = -p{D + FfD(pP/m n ) 
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~ Xp) /1 i \ 

(j (D + FfD(pP/m,) + -(D - 3F) 2 D(pf/m ,)) 

1?AA = gJip jjiD + ) 2 {v'fD(p { f /m K ) + pP D(pP/m K )) 

JP) 

hs+ = p( D ~ F) 2 D(p { p/m K ) 
hs°E 0 = _ F ) 2 (pp d (Pf ] / m K) + pf D(pf /mn K )) 

(n) 

>te- = - F) 2 D(pP/m K ) 

„ _ _ (-D — -F) (-D + 3F) / (n)p/ (n) / \ (p)t->/ H/ \\ 

r /SA - r/AE - Q7r 2p -P^IPf /"PC - Pf ^[Pf l m K )J 

where the function D(x) is 


( 21 ) 


( 22 ) 


For proton-neutron symmetric nuclear matter with Fermi levels pp = pp = 270MeV/c we 
hnd that 


AM* = AM; = — 154MeV , A M* K = -152MeV 
A M£+ = AMjo = AMj_ = — 18MeV 

A M| 0 = AMI. = 0 (23) 

and there is no 2° — A mixing. The values for AM* and AM* are in qualitative agreement 

with those extracted from nuclei (see, e.g. ||20|| ). For neutron matter with pp = 350MeV/c 
we hnd that the mass shifts are 

AM* = —57MeV , AM* = -lOOMeV 
A M* k = — 126MeV , A MZ+ = OMeV 
AMvo = — 15MeV , AM£_ = -30MeV 
AM| 0 = A Ml- = 0 . (24) 

The utility of our results depends crucially upon the smallness of higher orders in chiral 
perturbation theory at nuclear density. A computation at next order in chiral perturbation 
theory may help resolve this issue. We hope to present results on this in a future publication. 
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